function Tg = GGT(u)
global miu;
global IBd;

rI = u(1:3);
CbI = reshape(u(4:12),3,3);
rb = CbI*rI;

Tg = 3*miu/norm(rb)^5*cross(rb,IBd*rb);